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Abstract 

We review the critical behaviour of the Driven Lattice Gas (DLG) 
model. As a result, we obtain a novel Langevin equation for the DLG 
which depends on the microscopic transition probabilities. We then show 
how this dependence affects the critical behaviour of the the DLG plac- 
ing the finite and the infinite driving field cases into different universality 
classes. Two other well known anisotropic, conservative, non-equilibrium 
models, the two-temperature model and the randomly driven model (RDLG) 
are also studied. It is shown that the RDLG with infinite averaged field 
and the two-temperature model with Ty — oo fall in the same universality 
class as the infinitely driven DLG. A Langevin equation for the two-layer 
DLG is also presented. 

KEY WORDS: Non-equilibrium systems; Driven Lattice Gases; 
Langevin Equations 

1 Introduction 

Equilibrium statistical mechanics has succeeded in predicting the collective be- 
haviour of a system in thermal equilibrium given the laws governing its micro- 
scopic behaviour. Unfortunately, most natural phenomena belong to the field 
of non- equilibrium statistical mechanics, a subject far less well understood and 
still in its developing stage. Being overwhelmed by the enormous complexity 
of non-equilibrium systems, it is expedient to focus on those which settle into 
non-equilibrium steady states and to trace its behaviour to specific model in- 
gredients. Then, results pertaining to collective behaviour can be obtained by 
several methods. For instance, Monte Carlo simulations in which a great deal 
of our understanding of the large scale properties of a system is based and also 
serve to test predictions based on other approaches. Mean field theories, the 
usual starting point of the analytic route, can provide us with some insight into 



1 



the phase diagram, but they fail to describe the true collective behaviour prop- 
erly, specially in systems with local interactions in low dimensionalities. Fur- 
ther analytical developments are seriously hampered by, among other things, 
the (commonly) discrete structure of the system, so to describe the physics at 
the macroscopic level a new avenue is required. To accomplish this goal one can 
benefit from the mesoscopic approach in terms of Langevin equations, which 
concentrates on the long-time, large-distance properties and tries to eliminate 
the lattice altogether, by making the order parameter into a continuous field. 
This mesoscopic picture can be derived, at least in principle, starting from the 
microscopic master equation and implementing a coarse graining procedure, 
but in practice this route proves to be an insurmountable task for most sys- 
tems. This predicament is usually overcome by postulating phenomenological 
equations based on the choice of an order parameter and the underlying sys- 
tem symmetries. When equilibrium critical phenomena are considered, one can 
then appeal to the framework of the renormalization group, and universality 
then appears in the sense that the results are independent of the microscopic 
details, in particular the microscopic dynamic rules But non-equilibrium 
critical phenomena is still a challenging matter which displays striking features 
[^j. More specifically, in contrast to equilibrium, in non-equilibrium situations 
the transition rates are not a simply matter of convenience. The observable 
critical behaviour can depend on some details of the microscopic dynamics, a 
fact which is often underestimated in the literature. An analysis concluding 
about the relevant features that characterize the universal properties of a non- 
equilibrium system at criticallity is still lacking. 

This paper is devoted to the driven lattice gas (DLG henceforth) model and 
other related models. First devised by Katz, Lebowitz and Spohn [||, the DLG is 
a kinetic lattice gas of interacting particles subject to an external uniform field. 
A more detailed description of the DLG will be presented in the second section 
of this paper, but for the time being let us remark that various generalizations 
of the DLG have been argued to be relevant to the understanding of a wealth 
of varietes of natural phenomena (see |^ for a review). On the theoretical 
front, the DLG is one of the simplest non-equilibrium models that is hoped to 
serve as a paradigm for the behaviour of those systems which do not posses 
a thermodynamic equilibrium state. It also seems to capture the essence of 
strongly anisotropic systems. But before considering the DLG as a kind of non- 
equilibrium "metamodel" it is necessary to clarify some controversial issues. 
A longstanding glaring discrepancy exists on the value of the order parameter 
critical exponent (3. From Monte Carlo data for an infinite field (no jumps 
against the drive are allowed) a value close to 1/3 is found JH but, however, 
P = 1/2 is obtained from a mesoscopic Langevin approach ||^, ^ Several efforts 
aiming at solving this problem ||^, || have come to grief and the question remains 
under discussion. Besides this disagreement, the failure of the Cahn-Hilliard 
approach to coarsening dynamics |^ also casts some doubts on the reliability 
of the standard continuum approach to the DLG. Open questions abound and 
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clearly new inroads on this subject are needed. So, seeking a new framework for 
answers, we propose a new approach whose preliminary results were presented 
in0. ^ 

The remainder of this paper is organized as follows: In the first part of 
Section 2 we briefly summarize the basic ingredients of the DLG. We then 
give our derivation of a novel Fokker-Planck equation for the DLG and its 
stochastically equivalent Langevin equation. In Section 3 the two-temperature 
model, the Randomly Driven Lattice Gas and the multilayer variant of the DLG 
are studied employing the formalism of Section 2. Our concluding remarks are 
in Section 4. 

2 From Driven Lattice Gases to Driven Diffu- 
sive Systems 

Consider a set of particles confined to a box, A C Z'^ ^ with periodic boundary 
conditions. A configuration of this system is specified by giving all the site 
occupation variables, rix = 1, 0, echoing the fact that a particle may be present 
or not at site x. Besides this hard-core constraint, the model also includes a 
nearest-neighbour (NN) interaction so, given a configuration C = {nx}xeA, the 
Hamiltonian reads 

i?(C)--J^nxny. (1) 

NN 

The configurations C evolve according to a stochastic hopping dynamics which 
conserves the number of particles, or equivalently the density q. Up to the 
present, all we have is the familiar kinetic lattice gas for which the following 
master equation for the time evolution of the probability distribution Pt{C) 
applies 

dtPtiC) = {^[^' ^ C]PtiC') ~ W[C ^ C']Pt{C)y (2) 

c 

Here, W[C C] stands for the rate at which the system makes a transition 
from C to C", and C and C can differ by a single nearest-neighbour particle- 
hole exchange. The choice of W[C C] = D{PAH), where D is any function 
satisfying D{—x)= e^Z?(x) and AH = H{C') —H{C), ensures that the station- 
ary solution of (H) is the equilibrium one, i.e. P{C) oc e~^-^^^\ and /3 — 1/T is 
the inverse of the temperature of the thermal bath. 

Next, let us introduce a uniform (in both space and time) external drive E 
pointing along one of the principal axis of the lattice. We refer to it as the electric 
field while imagining that the particles behave as positive ions only in relation 
to it. The field biases the rates favouring jumps along its direction, suppressing 
jumps against it, and leaving unaffected those in the tranverse directions. For 
hard wall boundary conditions the only effect of the drive would be to add a 
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gravitational potential energy to the Hamiltonian (|I|) , and the resulting steady 
state would be an equilibrium one. But due to the periodic boundary conditions, 
the drive has a dramatical effect on the system static properties, preventing the 
system from achieving an equilibrium stationary state. In this case, the electric 
force is nonconservative so it is not derivable from a global potential. However, 
the local effect of the drive could be mirrored by adding to the Hamiltonian (Q) 
the work done by the field during the jump, so we choose the rates in the form 



W[C ^ C']= D{f3AH + f3eE), (3) 

where £ = (1,0, —1) for jumps (against, transverse to, along) E, and E is the 
strength of the electric field {E — |E|). When E = 0, one recovers the familiar 
equilibrium rates. 

Let us end this resume of the main features of the DLG (see |2[ |l^ for fairly 
detailed reviews) with a few words on its collective behaviour. At half-filling 
{g ~ 0.5) and _E = 0, a second order phase transition at the Onsager critical 
temperature, Tc, is known to occur in two dimensions. For E ^ and still at 
half-filled lattice, from the gleanings provided by Monte Carlo simulations, the 
DLG undergoes a second order phase transition at a higher critical temperature, 
Tc^\ saturating at about ri°°'' = lATc for E oo. The main observation is 

( E] 

that for temperatures above Tc the particles are distributed homogeneously 
while below Ti^^ the DLG segregates into a particle-poor phase and a particle- 
rich region, the latter having domain walls parallel to the field. 



2.1 The continuum limit 

In this paper we are mainly concerned with the continuum field description of 
the DLG, so firstly we aim at making the order parameter into a continuous 
field. To serve this purpose, let us mentally construct around each lattice point 
a boxi consisting of ri*^ sites. We then define a coarse-grained density (j)r,r G X 
through 

where f2r is the box centred at r. Now, we find inspiration in the original oc- 
cupation variables dynamics to postulate the time evolution of the new density 
variables: we perform exchanges between two randomly chosen nearest neigh- 
bour sites, r and r + a. Hence, if C""° is the configuration after the exchange, 
we have 

= {(/.X + r!-''(<5x,r+a - <5x,r)}xGA, (5) 

where the label a stands for the a direction. More generally, we can consider 
exchanges of magnitude rjQ.~^^ with probability amplitud f{rj), the latter being 
an even function of 77. Needless to say that 0r no longer equals 0,1, and when 17 
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is large enough the value of (pr becomes practically continuous, so that we have 

C''™ = {</>(x)+r;l]-'^V^„<5(x-r)KeT^ , (6) 

where T'^ is a d-dimensional torus. A suited Hamiltonian analogue to the inter- 
nal energy (pi) adopts the usual Landau- Ginzburg form 



H{C) = ^'^ j d'^x 



(7) 



The occurrence of the factor Qf^ in H{C) deserves some comments. It will 
play a decisive role when we discuss the influence of the transition rates in 
the dynamics. For the time being we shall just justify, albeit heuristically, its 
apperance. Let ^ be the system correlation length. Then, if we choose 1 ^ < 
^, the coarse-grained density will be aproximately constant over distances less 
than or equal to 17. After carrying partial summations over disjointed regions in 
A of linear size Vl, the factor Qf^ emerges and one can then take the continuum 
limit by simply making 17 large enough. 

Lastly, we associate a time dependent statistical weight with each configura- 
tion, Pt(C), which evolves in time accordingly to the following master equation 

dtPtic) = / (8) 

x|w^[C""''' ^ C]Pt(C"'™) - W[C C''™]Ft(C)}. 

The exchange rates W\C ^ C] have to be prescribed yet. We choose them 
such that 

W^p^C"™] =7:>(iJ(C"""'^)-i7(C)+il£;(C-^C"""")) , (9) 

with 

He{C ^C^n=r^8.-'E{\-4>{rf)+0{e) , (10) 

where e = 17"''. This means that the transition rates depend on the energy 
difference between configurations plus a term whose dominant part in e is the 
natural choice to mirror the effects of the drive as far as it accounts for the local 
increment of energy due to the driving field. Correcting terms of higher order 
in £ will be fixed later. Again, as in the DLG case described above, in absence 
of the drive the system stationary state is the equilibrium one characterized by 
the (f>^ Hamiltonian (^. 

2.2 Prom Fokker-Planck to Langevin 

We are part of the way towards deriving a Langevin equation for the DLG. 
The next step we take toward this end is to get a Fokker-Planck equation by 
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expanding the Master equation in e up to order. To keep matters simple, we 
will consider E=0 for the moment. We will recover the E ^ case later on. 

We avail ourselves of the results collected in appendix A to expand H{C^^°') 
and Pt{G'^'^°') around C. The Master equation transforms into 



dtPt{C) =J2j dvdinfiii) [{D{~^H) - D{AH))Pt{C) (11) 

5 

■ (50(r) 
where 

AH^ri\a-t^(v,^^K + 0{e^) , (12) 

with 

s:tt 

^a = -Vr^—- (13) 

We would like to stress that Aq is of order one rather than order e. This is 
due to the factor il'' in (|^) which ensues that our expansion of D is not around 
zero. This assurance is of the most importance because it bestows dependence 
upon the dynamics on equation (|l2|). Otherwise we would have found a very 
different story: our expansion would have resulted in a simple Model B [ p^ 
where any dependence on the dynamics would have vanished. Recalling again 
the results of Appendix A and noticing that much simplification is obtained 
because of the integration of odd terms in rj, one is led to 

dtPtiC)=J2 1 *(^-^) {ehiXa)PtiC) + '-^ei\,) (v.^J-^jPiC)}, 

(14) 

where 

hiXa) = f dv fiv) V DiV^a), 

Jr 

e{Xa) = [ df^ f{rj) rf D{t^K). (15) 

JR 

Turning our focus to the E ^ case, let us choose He in ( [lO| ) as 

i/.(C^C-) = ,A(f)+5:!L_J_£L(^Vr„Aj Ai-) (16) 



with aI^"* = a- E(l — 0(r)^). The point of these manouevres is that the Kramers 
Moyal expansion is now trivial. That is, one should only substitute in (llj) and 
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5|) Xa by Aq = Aa + aL^' . So, this election allows us to write 

5 



(17) 



+ye(A,) ( Vr 



(r) 



We have just derived a Fokker-Planck for the DLG. Now, we proceed to find 
out its stochastically equivalent Langevin equation. We shall invoke the main 
result of Appendix B in virtue of which the Fokker-Planck equation 



dtPtiC) 



{/,(</.;r)P,(C) + i (^Vr„^)5a(0;r)2Pt(C)}, (18) 



is equivalent to the Langevin equation 

d 



dt<j){r, t)^J2 /a(0; r) + ga{<l>; r)Ca(r, t) 



a=l 



(19) 



using the Ito prescription and with (a (r, t) being a gaussian white noise, i. e. (Ca (r, t)) 
and (Ca(r, i)Ca'(r',i')) = 5a,a'5{t — t')5{r — r'). Then, relating to our case, it 
is straightforward to get 



at0(r, t)=Y, Vr„ \h{k,) + e(A„)i/2^„(r, <) 



(20) 



where, time has been rescaled by a factor e, and finally e has been set to 1 since 
no more perturbative expansions in e are going to be considered. 

Before we proceed further, several comments are in order. First, the basic 
symmetries of the DLG are present in the Langevin equation (|^ : it is invariant 
under the simultaneous change E — > —E and — > — and it is also invariant 
under translations in space and time. But the central hallmark is that it depends 
strongly upon the dynamics. We believe this is a real step forward, if only 
because it goes beyond phenomenological approaches. Gratifyingly, we shall see 
shortly how to exploit this new state of affairs. 



2.3 Power Counting 

We focus on the critical region where large fluctuations on all length scales 
dominate. Further simplification in (20) is possible in this regimen by dropping 
the irrelevants terms in the renormalization group sense. Following the standard 
field theoretic methods let us introduce an external momentum scale /i and make 
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the following anisotropic scale transformations: i — > /i ^t, r± 



1 



and 



> H (j), where || stands for the direction parallel to the driving 
field E, and _L for those perpendicular to it. As usual, the noise scales as Ca 
^{z+d+cr~i)/2^^^ Next, we expand the Langevin equation in powers of fi around 
H = 0, keeping only the leading terms. The time scale, the transverse noise, 
and the transverse spatial interaction are forced to remain invariant under the 
transformation. With this understood the values of z and S can be determined. 
One gets z — i and 6 = {a + d — 3)/2. Different scenarios are now possible 
depending on the value of a. Demanding that the coefficients of and V| 
scale in the same way, as in the standard analysis of the DLG |^ , would lead 
to the choice a = 2. We show the most representative terms 



_e(0) 
2 

+h'{E) 



6" 



A 



-e(0) 



1/2 ' 



(3<T+rf-ll)/2£;y 



-ie(i?)i/2V||C| 



(21) 



which after setting cr = 2 and taking the limit /i — > 0, assuming d > 3, gives 



FTAi0 + fAi</)3 
6 



(22) 



h'{E) Vl(j}-Eh'{E) V||02 + ^V±C±(r) 



h' is the first derivative of the function /i(Aa), closely related to the first deriva- 
tive of the transition rate D. 

Let us take a glance at the structure of equation (22). It can be easily 



checked that all we have is the Langevin equation postulated by Leung and 
Cardy for the DLG and which is often known as driven diffusive system. But 
in stark contrast to equation (22) displays the precise form in which the 
microscopic field enters the mesoscopic picture of the DLG. More precisely, the 
two different critical temperatures introduced by Leung and Cardy for longitu- 
dinal and transverse ordering are indentified here as rh' [E) and r respectively, 
while the mesoscopic version of the field E finds its counterpart in Eh'{E). We 
summarize very briefly now the results of references ^: the upper critical 
dimension is shifted from the equilibrium value dc — A io dc = h. The reason 
for this change is claimed to lay in the strong anisotropic scaling jl^] (cr — 2) or, 
in other words, the necessity to introduce two length scales associated with the 
directions parallel and transverse to the field. An e-expansion can be therefore 
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performed in e = 5 — d dimensions — care must be taken of the quartic coupling 
g which is a dangerous irrelevant operator — leading to critical exponents that 
take their mean- field value (except the anisotropy exponent A This result 

is exact, i.e. higher order terms in the e-expansion vanish. This is borne out 
by a Galilean invariance symmetry present in the theory. In particular, the 
order parameter critical exponent /3 equals f/2. This figure does not tally with 
/3 w f/3 which is obtained from Monte Carlo simulations Q. Most of them, 
with few exceptions, use essentially infinite E, i.e. no jumps against the drive 
are allowed. It has been argued [|[ ^ that an anisotropic finite-size scaling 
analysis of computer simulations leads to critical exponents entirely consistent 
with field theoretic predictions. To date, the question remains murky so far an 
analysis of the data in |l^ was exhibited in [Q showing numerical evidence 
for /3 « 1/3. In the same vein, doubtful data collapse is obtained by all authors 
for the two point correlation function using anisotropic systems, specially at 
large distances Moreover, triangle domains pointing opposite to the field 
are found by numerically solving ( p^ ) while they are also found in microscopic 
simulations, but pointing along the field [^! Several models beyond the DLG 
have been proposed, practical motivations aside, to gain a wider perspective of 
these problems |l^ (see also later in this article). Despite these studies have 
led to new questions, some insight into the old ones have been achieved. That 
is, the very existence of a peculiar linear interface has emerged as a relevant 
point in the understanding of the DLG rather than the presence of a driving 
field. Also the existence of two different length scales is questionable p6[ |. 

It would be a mistake to dismiss all these puzzles which cast serious doubts 
on the validity of the existing field theory of the DLG. Originally, it was assumed 
that an infinite microscopic electric field implied a finite non-zero coarse grained 
driving field In equation (|2^) we show explicitely that this is not the case. 
In fact, it happens there that when E = oo the driving term disappears. We 
believe that this is the reason of the mismatch between the simulational results 
and the analysis of the Langevin equation (|2^). The latter is not suited for 
comparison with computer simulations because these are always performed with 
infinite drive and in equation (^) a finite non-zero driving term is present. Let 



us take the calculation a stage further by setting E to infinity in (22). Then, 
all the terms depending on the electric field E become identically zero. This 
fact can be easily checked, irrespective of the equation considered, (|2^) before 
power counting or (E2|) after rescaling. So, (E2h simplifies to 



Miv) = ie(0) 



-Ai</> + rA^<^+|A^03 
6 



(23) 



The latter equation constitutes a simple model B |13 in the transverse di- 
rections and no structure in the parallel one. Thus, the scaling a — 2 for E = oo 
leads to a trivial behaviour. Furthermore, one realizes that imposing and 
Vm to scale in the same way is meaningless because there is no parallel gradient 
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term appearing in the Langevin equation. Nevertheless, we have the freedom to 
choose a to look for different critical theories. A natural choice is cr = 1. With 
this election, an equation analogue to ( pT[ ) can be written down, however much 
more involved: 

_e(0) 
2 



9, 



h'{E) 
+h"{E) 



A±A|| 



^11 



.d-49 



|A||(/.3-^('^-8)/2i?V||(/.2 



/i('^-4)/V2V||(V|| 



)2 + /-41^£;A||03 + /i3('^-4)/2s2^||04 
o 



_^(5d-i6)/2^3:^^!Mv||,/,6 + e(0)i/2 y V^a + e(i?)^/'V||C|| . (24) 

± 

In particular, for E ^ one recovers the equilibrium theory (model B). For 
finite E, the critical dimension is dc = 8 but some of the irrelevant terms might 
be taken into account below the critical temperature. When E = oo we get a 
much more simple equation: 




-A±A|| 



Ai9i + TA_L(; 



^A, 



e(0) ^V^a(r,t) + V^ V||C||(r,i) 



(25) 



This equation is the central result of this work. We now proceed to discuss 
its physical implications. To begin with, equation ( p5| ) is structureless in the 
parallel direction. So, it corresponds to a lattice gas in which particles are 
exchanged at random in the direction of the field while jumps in the transverse 
directions are subject to energetics. The most obvious distinction between the 
finite E case and the infinite one lies in the steady-state current term. This 
does not appear in (^5|). Of course, such a current exists but it has no bearing 
on critical properties. Thus, equation (|25| ) does not gather it. 

It is quite surprising to find that, for infinite driving field, the Langevin 
equation changes dramatically when compared to the one that characterizes 
the finite E case. This would have been difficult to work out only on symmetry 
grounds. Also remarkable is the lack of Galilean invariance ||^ and the emergence 
of a single correlation length. Anticipating renormalization (the computation of 
the critical exponents will be presented elsewhere), we should remark that the 
upper critical dimension is now dc = 4 and it yields a distinct universality class 
from that obtained for finite fields. Therefore, a value of /3 different from 1/2 is 
expected. Finally, we wish to remark that equation ( p5| ) is renormalizablc when 
E = oo, while for finite values of the field it can be seen that it is not. A strong 
crossover from our theory to the finite E case occurs for very long impressed 
fields which obscures the interpretation of Monte Carlo data. 
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3 Related Models 



3.1 Two-Temperature model 

In this subsection we consider a model closely related to the DLG. The two- 
temperature Ising lattice gas consists of a set of particles endowed with an 
Ising Hamiltonian and lying on a hypercubical d-dimensional lattice. In con- 
trast to the usual kinetic lattice gas, particle-hole exchanges are controlled by 
rates with different temperatures according to the following rules: if the vec- 
tor a pointing from the particle to the hole lies in an n-dimensional subspace, 
the "parallel" subspace to say, then a transition rate of the form D(A_ff/T||) 
is at play. On the other hand if a belongs to the "transverse" space then the 
exchanges are coupled to a bath with temperature T±. Although this model 
could be seen somewhat afield from the DLG — there is no driving field, for 
example — Monte Carlo data for both the DLG and the two-temperature model 
for T|| = cxo look quite similar jl^. For instance, the value for the order parame- 
ter critical exponent in the two-temperature model is /? « 1/3, thus sharing the 
same value as the DLG with E = oo. To illustrate this connection further, let 
us consider our continuum version of the DLG with infinite drift. As we saw in 
the previous section this corresponds to completely random particle hops in the 
field direction. Then we note that the same picture arises if we choose T|| = oo. 
So, we should expect that the two models lie in the same universality class. 

We now proceed to build a continuum theory for the two-temperature model. 
Let us start with the master equation in continuous space (||). Now W[C ^ 
C] = D{H{C') - H{C)) where H{C) is the usual 0-four hamiltonian (|) but it 
depends on a that the coefficient t of (/'^/2 will adopt different values. That is, 
r = if a lies in the transverse subspace and t = ry otherwise. Being initially 
interested in the t\\ — > co case, we shall assume that particle-hole exchanges in 
the parallel subspace are subject to 

H ^n'^^- j dx02(x). (26) 

So, it is expedient to separate (^) into two parts, the first being a sum over the 
transverse subspace and the second over the parallel one. Next, we follow the 
same steps as in Section 2.2, i.e. an expansion of AiJ and Pt{C''''^°') around C. 
The calculation is straigtforward, for the (simboUicaly) a G-L case leads to the 
same result as in Section 2.2 and when a e|| one simply gets A_ff = — 7/r|[ Vr„<?!>. 
Then the following Fokker-Planck equation can be easily computed 

BMC) = Y^Jdv{^^,,j^^^\ehiX)P,{C) 
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aG II V 



(r) 



PtiC), 



(27) 



and in the limiting case r|| — > oo its equivalent Langevin equation reads 



(28) 



as 



Finally, we restrict ourselves to a 1-dimcnsional parallel subspace. A naive 
dimensional analysis yields, after dropping irrelevant terms and assuming a = 1, 
the same equation we found for the DLG at criticallity with E — oo. Therefore 
we are supplying a firm grasp to support our claimings. That is, the DLG with 
E = oo and the two-temperature model with T|| — oo are members of the same 
universality class. 

No adiditional effort is necessary to study the finite Tn case. We again seper- 
ate (^) into two parts, but now we take the full expression for the Hamiltonian 
(^ rather than (p6|). A derivation running along the same lines of the previ- 
ous subsection is possible. We only remind that the scaling z — 4,(7 — 1 and 
S = (d — 2)/2 is used in the power counting procedure which entails 



e(0) 
2 



A20+(riAi+T||A||),/)+|^A(/.3l +y^^Vr„Ca(r,i)- (29) 



Despite of the fact that anisotropy is only present as far as the masses are 
concerned, we should remark that more general anisotropics could be expected. 
Then, the collective behaviour of the two-temperature model with finite T|| could 
be more adequately predicted by means of an extension of ( p9|) to an equation 
with full anisotropy in the coefficients. 



3.2 Randomly Driven Lattice Gases 

Here we take up an extension of the DLG, the randomly driven lattice gas 
(RDLG hereafter) . Let us consider a DLG in which the driving field fluctu- 
ates accordingly to an even distribution p[i?(x, t)] which is (5-correlated in space 
and time. The RDLG is easier to realize in the laboratory than the DLG and 
one can also benefit from a higher analytical simplicity because a random drive 
induces no steady-state current and the particle-hole symmetry is preserved. 
Like the DLG, it exhibits a second order phase transition at half filling from 
a disordered state to striplike order. There is hardly any difference between 
typical ordered configurations associated with the DLG with an infinite drive 
imposed, the two-temperature model and the randomly driven model. Now, to 
provide a better comparison we construct a continuum equation for the RDLG. 
Our starting point is again equation (||) with the prescriptions (|^) and (|lO|) . We 
should then average over a random E, but we can defer such an average to a 
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later stage. First, we repeat the steps we performed to arrive at a Langevin 
equation for the DLG. All results carry over without change until we face to 
averges over E in ( p^ of the type E"^h''"'\E) and y/e{E). As a consequence 
of a symmetric ^[-^(x, t)] and the integrations over 77, the averages with m + n 
being an even integer vanish whilst the remainder terms will yield the following 
finite values: 

71 = J dE p[E]h\E), 

J dE p[E] E h"{E), 

J dE p[E] y/e{E). (30) 



72 
73 



The Langevin equation can then be written down in the form 

dtcj, = -^5] Vl0 + 7iAf + (71 - fM)A^A(|0 + ^rAi0 - r7iA„<^ 

+ ^|a^(/.3 + ^r(72 -7i)A||03 ^ ^^^^^^^ (3^) 

We have arrived at an entirely anisotropic equation, i.e. all the gradient opera- 
tors have been split into components parallel and transverse to E. In accordance 
with our earlier discussion of the two-temperature model, we conclude that the 
RDLG and the two-temperature model with finite T\\ share the same critical 
behaviour. But we caution that, as in the DLG, in the RDLG two cases have 
to be distinguished. We can think of the simplest distribution p[E], namely the 
bimodal ^[5{E-\-Eo) + 5{E — Eq)\, and then take the limit Eq — > 00. Not surpris- 
ingly, in this limit 7i = 0, i = 1, 2, 3, and the Langevin equation ( |25| ) associated 
with the infinitely driven DLG emerges again. We therefore can conclude that 
the generical critical properties of the two-temperature model with T|| =00, the 
DLG with E set to infinity and the RDLG driven with an effectively infinite 
field are indistinguishable. 



3.3 Layered Driven Lattices Gases 

We turn our attention to a generalization of the DLG Let us consider a 

pair of identical square lattices placed back to hack. Each plane is a copy of 
a two dimensional DLG. No inter-layer coupling is allowed, but particles can 
hop from one plane to their "nearest neighbour" site in the other one. This 
process is controlled by in-plane energetics alone and it is not affected by the 
drive. The overall particle density is fixed at 1/2. Concerning the nature of 
the phase transition, Monte Carlo data [E = 00) have revealed that intriguinly 
two transitions appear j2^. Here, we just give the gist of references |7[ ^ 
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where the phase diagram in the {E, T)— plane has been mapped out employing 
Monte Carlo simulations and dynamic mean-field theory. As T is lowered the 
system first orders from a homogeneous state into a state with strips in both 
layers. This transition is observed to be characterized by the same critical 
indexes as the DLG. Decreasing T further, we reach the second transition where 
homogeneous layers with different density appear. This transition belongs to 
the Ising universality class for any E < E^. ^ 2, while values of E beyond the 
threshold field E,, lead to a first order phase transition . 

Although physical motivations for this model come from various directions 
we shall focus on the theoretical side. As far as comparison with simula- 
tions is intended, we shall provide a mesoscopic picture to place into a coherent 
analytical context these two phase transitions. Next, we construct a continuum 
mesoscopic theory in the spirit of Section 2. We first discuss the equilibrium 
case E = 0. Jumps in this model can be naturally divided into two types: 
m-layer jumps and inter-layer jumps. So, the following notation will prove to 
be convenient: we shall refer to 4>i{^) and 02 (r) as the coarse grained density 
field in plane one and two respectively, and an arbitrary global configuration 
will be termed C = {(j>i (r), (/)2(r)}. We shall denote C^'^"' the configuration after 
an exchange of density erj is performed in the a direction with an infinitesimal 
neighbour of r in plane i. Exchanges between planes will lead to configurations 
named C"". More specifically, 

Cr = {0i(x)+£77Vx„<5(x-r),02(x)}, 

Cr = {0i(x),02(x)+£77Vx„<5(x-r)}, 

C"- = {0i(x)+£775(x-r),(/.2(x)-£77<5(x-r)}. (32) 

With this understood, then the following master equation can be written down, 

dtPtiC) = E /^^ dr^f{v){w[Cr ^ C]Pt{Cr) " W[C ^ Cr]Pt{C) 

+W[CT" ^ C]Ft(Cf ") - W[C cr"]-Pi(c)} 

+ J drdrj f{rj) |p^[C"'''" ^ C]Pt(C"'''") - W[C -> C"'''"]Pt(C)}. 

(33) 

W[C C] has its usual meaning, i.e. W[C C] = d[h{C') - H{C)^ 

and H{C) is again the Hamiltonian (^. The calculus towards a Fokker-Planck 
equation can be carried out as we did in Section 2, the only difference being the 
inter-layer current term. One can easily get 

dtPiC) -t.T.Jdr (^-^) {eMA«)P*(C) (34) 
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We now explain our notation. The first term in this equation models two de- 
coupled Ising lattice gases with appropiate constraints. As for the second one, 
it is simply a "discrete" version of the former, as long as it has to be with ex- 
changes across the layers. Thus, Xa has the same meaning as in Section 2, but 
restricted to plane i. The functions h and e are also defined as in the previous 
section. With Vi2^ we simply denote the operator 



and Ai2 stands for ^Vi2^^ H{C). Turning the drive on, one only has to move 

Aa'* into Aq-* thereby taking into account the effect of the electric field. A 
Langevin equation can then be derived following the same lines of Subsection 
2.2. One gets 

^^4>^{v) = -rAi2 - re(Ai2)^/'C(r, t) + J] V^, [h{K(^'>) + e{K(^yf'\,{v, t) 

a 

5t02(r) = rAi2+re(Ai2)i/2C(r,i)+5]Vr„[/i(Ai2)) + e(Ai2))V2^4r,t) 

a 

(35) 

^ and Ca are gaussian white noises while F is a hand introduced transport 
coefficient that measures the rate at which the system changes due to the inter- 
layer exchange mechanism. 

Bearing density conservation in mind, we introduce two new fields: 

77l(x) = ((/.I + (/.2)/2, </j(x) = ((^1 - (/)2)/2. (36) 

Equations ( ^5|) are easily expressed in terms of the new fields m(x) and (^(x). 
Now, simulation results hint at which field will be treated as an order parameter. 
We shall take m(x) as the ordering field for the DLG type transition, the one that 
occurs at a higher temperature. Then, in a naive dimensional analysis there is 
much freedom to choose the scale of observation. We perform the following scale 
transformation: t —> fi~^^t,r — > /i^'^r, — > and m — > ji^m. In particular, 

if we fix the exponents z — 4,a = 1,S ~ d/2 and 7 = (d — 2)/2, it can be easily 
checked that, after neglecting terms that are irrelevant in the renormalization 
group sense, we are left with nothing but equation (^) for the DLG. Hence, 
the critical properties belong to the DLG universality class, a picture consistent 
with simulations. 
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Turning next to the second transition, we consider tp as a non-conserved 
order parameter. We propose a critical theory naively consistent with z = 
2,cr = 1,<5 = (d — 2)/2 and 7 — d/2. This scaling leaves us with a couple of 
equations that take the form: 

dtm = r(/i'(0)A_L + /i'(£;)A||)m-/-2|(/i'(0)Ai + ;i'(i;)A||)(^2^) 

+ VW)Y1 ^^(Cl'^ + Cf ) + V^V|| (C,;') + C|f (37) 

We have dropped all terms that give a negligible contribution in the limit /i — > 0. 
Due to the electric field, which singles out a lattice axis, all gradient terms have 
become anisotropic. The situation is then very much reminiscent of the driven 
lattice gases with repulsive interactions We have an electric field E that 

has an effect on the phase transition only through an auxiliary non-ordering 
field to(x). The naive dimension of E turns out to be (d — 2)/2, in contrast 
to (d — 8)/2 (Section 2.3), so it is highly irrelevant compared to gtp^ . We note 
that essentially the same set of equations results in |2^ , so our analysis of ( ^7|) 
will follow the same lines of this reference. That is, E is naively irrelevant for 
the Gaussian fixed point until d is lower than two. We conjecture that E is 
not relevant to the Wilson-Fisher fixed point for d < 4, so the unique effect of 
the drive consists in generating anisotropics. It should be noted that the field 
m(x) does not order, so we do not need to keep track of it as far as critical 
behaviour is concerned. Thus, the critical properties for the two-layer driven 
lattice gas are given by the Langevin equation for iy9(x) and they fall into the 
Ising universality class. However, corrections to order 0{E^) show that g and t 
decrease to an amount that depends on E. Eventually, both of them may vanish 
simultaneously, a mechanism that would be liable for a tricritical point. Then, 
in qualitative agreement with simulations, the transition would be discontinous 
for E beyond a critical field value. Interestingly enough, it could be worked out 
the dependence of the transition temperature on the dynamics. 

4 Summary and Conclusions 

Here, we take stock of what we have done. Section 1 outlines the task under- 
taken in this article: to arrive at a field theoretic description of the DLG model 
and to show how this description explains why the DLG, and three models in- 
timately related to it, have the various properties that they exhibit. Inspired 
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by the dynamics at the microscopic level, we were able to construct a master 
equation in continuous space. The goal of the avenue we have pursued was to 
take into account the microscopic details, an effort towards disentangling the 
role of dynamics in non-equilibrium critical behaviour. An expansion of the 
continuous master equation was then possible due to the factor r2'^, a remnant 
of the coarse graining over the underlying lattice. The role of the factor can 
only be fully appreciated noticing that all dependence on the dynamics relies 
on it. Then, the full Langevin equation for the DLG was derived. The richness 
of equation ( po| ) was brought to fruition when viewed on different length scales, 
leading to the emergence of distinct critical theories. In particular, we have 
recovered the Langevin equation of reference j|] after considering an anisotropic 
scale transformation (cr = 2). This result, interesting though it is, is unsatisfac- 
tory in several respects. Firstly, simulations favour a value /3 1/3 for the order 
parameter critical exponent, whilst a mean field value /3 = 1/2 ensues exactly in 
for 2 < d < 5. Diniing further the validity of equation (|2^) as a continuum 
description of the DLG, numerical work on equation (|2^) conflicts with the sim- 
ulation data The situation became more transparent after we resorted to 
the detailed dependence of the coefficients on the microscopic dynamics. That 
is, the case E = oo, which is with few exceptions the most studied one, was 
carried out explicitely the result being a trivial equation. This is a consequence 
of the choice cr = 2, which makes no sense as it was seen in Section 2. Following 
well honed arguments ||l^ that invoke a single effective correlation length, we 
have turned to the choice ct = 1. Possibly against intuition, different critical 
behaviour has been found for the finite E and E = oo cases. The upper critical 
dimension associated with the former is dc = 8 whereas the latter is character- 
ized by dc = 4, thereby in either case yielding an universality class other than 
that obtained in 

In a subsequent section we have provided three examples of the applicability 
of our methods in the shape of the two-temperature model, the random driven 
lattice gas and the driven bi-layer lattice gas. On symmetry grounds we have 
suggested why the DLG and the two-temperature model with T\\ — oo should 
lay in the same universality class. Our suspicions have been confirmed so far a 
Langevin equation identical to the one associated with the DLG has resulted in a 
derivation that runs along the same lines of Section 2. The finite case has also 
been studied. Different critical behaviour from the T|| — oo case has been found. 
Remarkably, the RDLG model with finite averaged external field has resulted 
in a Langevin equation identical to that associated with the two-temperature 
model with T|| finite, whilst the RDLG with an infinite averaged driving field 
was proved to belong to the same universality class of the infinitely driven DLG. 
Summing up, we conclude that the infinitely driven DLG, the two-temperature 
model with T|| = oo and the RDLG with infinite drive, are described by the same 
Langevin equation. On the other hand, the two-temperature model with finite 
T|| and the RDLG with finite averaged field belong into the same universality 
class. 
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In the last part of Section 3 we have tackled with the two-layer driven lattice 
gas. A viable explanation for the two transitions exhibited in this model has 
been provided in the frame of field theory. Again, the election a = I has proved 
to suffice our purposes, i.e. the understanding of collective behaviour in these 
systems. Clearly much work is still to be done as more detailed studies on these 
subjects would be highly desirable. 



A Formal developements 

Supose that F is a functional of a function (j){x). If (p changes to + 50, then 
a Taylor series expansion can formally be written down 



F{cl> + 5ct>) =F{<i>) + [dx, JI—5ct>{x^ ) 

J 0(p[Xi) 



(38) 



where 5F/5(j) means the functional derivative of F{(f)) with respect to (l){x). Now, 
we explicitely treat the case we are concerned with, namely 

5(i) = EVJ{x-r). (39) 

In such a case it inmediately follows that 

F{<t> + 54>) = F{<t^) + f]l- f dx,...dxn ^r^^tL A 
^n!i 6(j}{xi) . . . d(j){xn) 

X£"Va;i(5(a;i -r)...VxJ{xn -r) 

oo „ J- n „ A 'I 



Tl=l 

The operator (V^^) satisfies 



\ 5(j){r) J 5(j}{r) 

which can be better proved by putting the last expression in a lattice. Finally, 
the usual properties of functional derivatives can be applied, v.g. 

V.A)^,(^)^,(^)=^,(^)(v.^)f.(0) + F.(0)(v.^)f,(0). 

(42) 
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B Fokker-Planck and Langevin equations for sys- 
tems with conserved order parameter 

A Langevin equation with conserved order parameter has the general form 

dt^{x,t) =J2^^a [fa{^,x)+ga{^,x)Ux,t)] , (43) 

a 

(a being a gaussian white noise: {C,a{x,t)z'g^{x' ,t')) = 5a,a'5{x — x')5{t — t') ; 
{C,a{x,t)) = 0. Let us introduce a time discretization. Then 



^+1 = (l>n + [fa{(l>n,x) + ga{<Pn, x)Ca,: 



(44) 



where (/>„ (f)t and tn = ne ^ t when e ^ and n — > oo. A factor e has been 
absorbed into the noise and we are using the Ito prescription. 

The probabihty that the system is in the configuration (/)„+i(r) at time tn+i 
is given by 

P„+i(</.„+i) = (y #„P„(</.„) S{Sn))^, (45) 



S{Sn) ~ 5{(f)n+l -(f)n) + j dxi 

if 



66 {(pn+l - 4'n) 



+ 



6H{(j>„ 



+ - / dXidX2 . , , , 
2 7 64)n\X\)64)n\X2) 



a a 

The calculus now reduces to noise averages. An intermediate step is 

\PnY,^ 6{x^-X2)) 



+ ... (46) 



- / dx\dx2 



6(j){xi)6(t){x2) 



(47) 



and after a bit of algebra, applying the results of Appendix A, our final result 
reads 



w^)=5:/dx(v.„^)[pj„+i(v^ 



' 6(pix) 



(Pt gl) 



(48) 
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